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ABSTRACT 

We present here a direct comparison of the polarisation position angle (PA) profiles of 
17 pulsars, observed at 1.4 and 3.1 GHz. Absolute PAs are obtained at each frequency, 
permitting a measurement of the difference in the profiles. By doing this, we obtain 
more precise rotation measure (RM) values for some of the pulsars in the current 
catalogue. We find that, apart from RM corrections, there are small, pulse longitude 
dependent differences in PA with frequency. Such differences go beyond the interpre- 
tation of a geometrical origin. We describe in detail the PA evolution between the two 
frequencies and discuss possible causes, such as orthogonal and non-orthogonal polar- 
isation modes of emission. We also use the PA and total power profiles to estimate 
the difference in emission height at which the two frequencies originate. In our data 
sample, there are changes in the relative strengths of different pulse components, espe- 
cially overlapping linearly polarised components, which coincide with intrinsic changes 
of the PA profile, resulting in interesting PA differences between the two frequencies. 

Key words: pulsars: general - polarisation 



1 INTRODUCTION 



discovery of pulsars, 
developed a model of 



Shortly after the o ri ginal 
iRadhakrishnan fe Cooka |l96£T 
emission which has subsequently been used to interpret pul- 
sar polarisation. According to this model, pulsar radiation 
is polarised along, or orthogonal to, the open magnetic field 
lines of the strong dipolar magnetic field that surrounds 
the star. As the beam of the pulsar sweeps past the line of 
sight of the observer, the position angle (PA) of the linear 
polarisation rotates, resulting in the well known S-shaped 
curves of PA versus pulse phase seen in many pulsars. 
The angles which determine the sweep of the PA are the 
inclination angle a between the magnetic and rotational 
axes and the impact parameter /3, which is the smallest 
angle between the locus of the line of sight on the pulsar 
beam and the magnetic axis. This rotating vector model 
(RVM) has been used in some cases to infer the geometrical 
angles of pulsars from polarisation observations. 

Despite the initial success of the RVM, numerous at- 
tempts to fit polarimetric data demonstrated that in many 
cases, the parameters of the fit were poorly constrained, 
mainly due to the restricted longitud e over which pulsar s 
emit. For instance, in the case of lBlaskiewicz et alJ lll99ll) . 
who are primarily concerned with the point of steepest gra- 
dient of the PA swing, the fitted parameters a and /3 at 430 
and 1418 MHz rarely agree. The authors attribute this to 



very large error estimates, often as large as hundreds of de- 
grees. The same is seen in anot her set of publ ished RVM fits 
at mu lt iple radio frequen cies llvon Hoensbroech fc Xilourisl 
ll997ah . lMitra fc Lil <2004l) recently presented RVM fits for 
6 pulsars, to data obtained at multiple f requencies (4 to 6 
freque ncies for each pulsar). Similarly to Blaski ewicz et ail 
( 1991), the fitted parameters (especially a) generally do not 
agree between the different frequencies, despite the claim 
of the aut hors to the contrary . In a remarkably complete 
treatment. lEverett fc Weisberd (1200 if) summarise the prob- 
lem and demonstrate the difficulties in obtaining a unique 
RVM solution. They attempt RVM fits on their data, ob- 
taining reliable results for only 10 out of the 70 pulsars in 
their sample. Even in the limited cases, the fitted parame- 
ters are not always in agreement with previously published 
results. 

Early observations also showed that average pulse pro- 
fil es become narrower with increasi ng frequency. The model 
of iRuderman fc Sutherland! (Il975l) suggests that the fre- 
quency of emission is related to the local plasma density, 
which drops off at higher altitudes from the surface of the 
star. This leads to a natural radius-to-frequency mapping 
(RFM), with higher frequencies originating closer to the 
pulsar surface. In the dipolar magnetic field structure, the 
open field lines spread out at greater altitudes, ex plaining 
the broader widths of profiles at lower frequencies llCordej 
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Il978l) . Although the exact height at which a particular fre- 
quency is e mitted depend s on the nature of the emission 
mechanism (|Melrosei 12000). some sort of RFM appears to 
underly the observations. 

To first order, the PA at a given pulse longitude is 
fixed, irrespective of the emission height. The consequences 
of RFM on the PA lie mainly in second order effects, 
such as relativistic effec ts within the pulsar light cylinder. 
Blas kiewicz et al] Jl99lT) showed that due to such effects, the 
centroid of the PA swing is delayed with respect to the cen- 
troid of the total intensity by an amount proportional to the 
radius at which the radiation originates, providing a method 
to estimate relative emission heights in the magnetosphere. 
Such calculations agreed well with the heights inferred by 
the changing widths of average profiles. 

In a c omparison of data from PSR B0329+54 at two 
frequencies lMalov fc Suleimanoval (1199a) showed that the 
position of the pair of outer components changes relative 
to the central component, due to effects of retardation 
and aberration. A method for estimating emission altitude s 
was subsequently proposed bv lGangadhara fc GuptalteOOl]) . 
who exploit the fact that components originating close to the 
magnetic axis are emitted from much lower altitudes than 
pairs of components that flank them. In pulsar profiles with 
an appropriate configuration, the relativistic phase shift be- 
tween the central (core) component and the middle of the 
outer component pair provides the only information nec- 
essary to compute emission hei ghts, as demonstrated in a 
refi nement of the me thod of iGangadhara fc Guptal (1200 1]) 
by iDvks et al] (I2004) . The emission altitudes using this 
relativistic method are generally larger than those of the 
Blas kiewicz et al] (ll99lT) method. In practice, the main dif- 
ficulty of estimating emission heights in this way relates to 
the limited pulsar profiles with a central component and 
outer components that can be unambiguously identified as 
a pair. 

There are several reasons that cause observed PAs 
to diverge from the RVM. The first and probably best 
documented phenomenon is ortho gonally polarised modes 
(OPM) in pulsar emission (e.g. [ Manchester et al] ll9T5t 
IBacker et alll975t ICordes et alJl97sl) . These result in parts 
of the PA profile being 90° offset to the main swing, an 
effect which can be identified and accounted for. Also , 
OP Ms are frequency dep endent fe.g. lStinebring et al.lll984l 
iKarastergiou et al]l2002f) . which results in different PA pro- 
files at different frequencies, the higher frequencies showing 
more occurrences of 90° discontinuities. 

An additional source of the PA dependence on 
frequency may originate from non-orthogonal polarisa - 
tion modes, identified first by IBacker fc Rankira il98Cf) . 
iRamachandran et al] i2004h recently showed that the dif- 
ferent spectral behaviour of quasi-orthogonal modes in PSR 
B2016+28 results in a ch ange in the PA swing over a very 
narrow frequency range. IKarastergiou et al] i2005h claim 
that the spectral behaviour of the total power and linear 
polarisation in different components can also be attributed 
to different spectra of the orthogonal modes. The possibility 
of variable PA profiles at different frequencies then naturally 
arises. 

iMichell <ll99l]) considered Faraday rotation in the mag- 
netospheres of pulsars as the cause of complexity in the 
PA profiles of pulsars. However, he argued that the phys- 



ical grounds for this process are lacking and subsequently 
no attempt was made to observationally identify this phe- 
nomenon. 

All of the above mechanisms that distort the simple 
geometric nature of the RVM, also have some form of in- 
herent dependence on frequency. We therefore proceeded to 
investigate this by presenting for the first time direct com- 
parisons of the PA swings of 17 pulsars between 1.375 and 
3.1 GHz. In doing so, we document the nature of the fre- 
quency evolution of the PA, which we discuss in the context 
of the above physical processes and models. The purpose of 
our analysis is to specifically avoid RVM fits to our data 
and geometrical interpretations which rely upon these fits. 
On the contrary, instead of comparing the PA profiles to a 
badly constrained curve, we compare them as precisely as 
possible to each other in an unbiased approach to study PA 
changes. We believe this approach is warranted given the 
nature of most PA profiles in our data and the difficulties in 
fitting RVM curves over very small longitude ranges. 

In Section 2 we describe the observations. In Section 3 
we present a method by which we align the profiles at the 
two different frequencies using the absolute PA values, which 
also permit us to determine very precise RMs. In Section 4 
we discuss the polarisation profiles and compare them to 
previous profiles at other frequencies. In the final parts of 
this paper, we discuss the differences between the profiles 
at the two frequencies and possible interpretations in the 
context of the above models for the origin and frequency 
evolution of PA profiles. 



2 THE OBSERVATIONS 

The observations were carried out using the Parkes radio 
telescope on 2004 August 31 and September 1. We used the 
H-OH receiver at a central frequency of 1.375 GHz with 
a bandwidth of 256 MHz and the 10/50 cm receiver at a 
central frequency of 3.1 GHz with a bandwidth of 512 MHz. 
In both cases, the backend correlator subdivided the total 
bandwidth into 1024 frequency channels and also recorded 
1024 longitude bins per pulse period. The receivers have 
orthogonal linear feeds and also have a pulsed cal signal 
which is injected at a position angle of 45° to the two feed 
probes. 

The pulsars were observed for 30 minutes at each fre- 
quency. Prior to the pulsar observation a 3 min observation 
of the pulsed cal was made. The data were written to disk 
in FITS format for subsequent off-line analysis. Data anal- 
ysis was carried ou t using the PSRCHIVE software package 
jHotan et aTll2004h . 

We used standard techniques to calibrate the relative 
gains and phases of the two linear feed probes. The pulsed 
cal is injected at 45° to both probes so that its intrinsic 
signal can be assumed to be 100 per cent polarised in Stokes 
U. Gain corrections can be made by comparing the observed 
Stokes Q with Stokes / and phase corrections can be made 
by comparing the observed Stokes U with Stokes V on a 
channel by channel basis in the output data. Following this 
calibration Q and U can be further rotated to correct for the 
angle which the feed probes make to true north and for the 
parallactic angle of the source at the time of the observation. 
At the end of this procedure we have the intrinsic PA of 
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Name 



Period 



W10 



RM 



RM reference 



J2000 


B1950 


111s 


1/1 


' 7 2 


rad/ 










(Ic <->' 


d £' 


New 


Previous 


J0738-4042 


B0736-40 


375 


32 


30 


14±0.5 


14 


J0742-2822 


B0740-28 


167 


16 


16 


150±0.5 


156 


J0835-4510 


B0833-45 


89 


14 


16 


32±0.5 


38 


J0837-4135 


B0835-41 


752 


5 


6 


146±1.0 


136 


J0942-5552 


B0940-55 


664 


26 


29 


-62±0.8 


-62 


J1056-6258 


B 1054-62 


422 


32 


32 


-1 or 6±0.5 


4 


J1157-6224 


B1154-62 


401 


45 


48 


507±2.5 


508 


J1243-6423 


B1240-64 


388 


7 


8 


162±1.0 


158 


J1326-5859 


B1323-58 


478 


17 


16 


-586±0.7 


-580 


J1327-6222 


B1323-62 


530 


11 


11 


-319±0.8 




J1356-6230 


B1353-62 


456 


32 


27 


-580±1.7 




J1359-6038 


B1356-60 


127 


14 


14 


39 or 37±0.5 


33 


J1602-5100 


B1558-50 


864 


11 


13 


79±1.3 


72 


J1644-4559 


B1641-45 


455 


14 


11 


-618±0.5 


-611 


J1752-2806 


B1749-28 


562 


8 


8 


91±3.0 


96 


J1807-0847 


B1804-08 


164 


28 


28 


171±1.8 


166 


J1825-0935 


B1822-09 


769 


22 


21 


68±0.7 


65 



Taylor et al. (1993) 
Han et al. ( 19991 
Hamilton et al. (1977) 
^^ylor^^aL 
Tavlor^^aL 
Costae^al. 
Tavlor^^al. 
Tavlor^^aL 



Costa^^ah 



(1993) 
' (1993) 
(1991) 
(1993) 
(1993) 
(1991) 



Han et al. (1999) 
Tavlor et al. (1993) 
Taylor et al. (1993) 
Hamilton k, Lvne (1987) 
Hamilton fc Lvne (1987) 
Hamilton fc Lvne (1987) 



Table 1. The list of pulsars and their parameters. v\ and V2 correspond to 1.375 and 3.1 GHz respectively. 



the radiation just above the feed at a particular frequency. 
PA values are only computed when the linear polarisation 
exceeds the noise by a factor of 5. 

To obtain the PA at the pulsar, the Faraday rotation 
both through the Earth's ionosphere and the interstellar 
medium must be corrected for. This is done in a two stage 
process. First, we fit for the RM across the 256 MHz band 
at 1369 MHz. This yields an RM with an error of order 
±5 rad m~ 2 . The PAs at both 1.375 and 3.1 GHz can be 
rotated by an amount given by 

(1) 

where c is the speed of light and v the observing frequency. 
A more precise RM can then be produced by comparing 
the PAs at both frequencies, using a least squares fitting 
technique as described in more detail below. 



where ctpai and erpA2 are the uncertainties in the PA mea- 
surements. As shown later, the instances of orthogonal PAs 
at the two frequencies are few and do not influence the value 
of Xdx ■ Despite this, if any instantaneous difference in the PA 
is greater than 45° , we assign the PA difference the comple- 
mentary value to 90° , thus avoiding great dPA values in the 
case of orthogonal PAs at the two frequencies. We then shift 
the profile by one bin and repeat the process, thus yielding a 
table of x 2 as a function of the shift dx. The value of dx that 
corresponds to Xmin ls chosen as the shift in pulse-longitude 
to align the two profiles, and the corresponding average dPA 
is used to apply a vertical shift. 

The vertical shift that is required to align the profiles 
can be achieved with a correction in the RM value and there- 
fore by calculating the mean dPA offset we are computing a 
very accurate RM. The error on the RM then corresponds 
to the error in the mean dPA, which is 



3 DATA ANALYSIS 

A new approach is used in this paper to align profiles at 
1.375 and 3.1 GHz. Initially, we slide one PA profile over the 
other in the horizontal, pulse-longitude direction, in steps of 
one bin of data. At every step of shift dx, we measure the 
average offset dPA between the two PA profiles 



OdPA 



-Y 



(4) 



The contributions to the measured values of RM aris- 
ing in the ionosphere have been estimated by integrating 
a time-dependent model of the ionospheric electron density 
and geomagnetic field through the ionosphere along the sight 
line between the telescope and the pulsar. The average con- 
tribution is found to be ~ — 1 rad m -2 . 



dPA d 



IV 

l^PAl i -PA2 i 



(2) 



where N is the overlapping number of data bins which have 
a PA measurement at both frequencies, i the index of each 
such bin and PA1 and PA2 the value of the PA at that bin 
at 1.375 and 3.1 GHz respectively. Then, we calculate a nor- 
malised chi-square between the two sets of data, considering 
their average PA offset, according to the equation 



1 ^ (PAli -PA2i -dPA dx ) 2 



Xdx / J 



2(a| 



(3) 



4 RESULTS 

Table 1 shows the list of sources, their rotation period, the 
pulse widths at the 10% level of the peak flux density at both 
frequencies and the RM, after removing the ionospheric con- 
tribution. For comparison, we provide the current published 
RM values where available. 

For an accurate comparison between the polarimetric 
profiles at 1.375 and 3.1 GHz, we have overlaid them for 
each pulsar, as shown in Figure 1. Due to the steep flux 
density spectrum of pulsars, the sources are much brighter 
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PSR J0738-4042 
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Pulse Longitude (deg) 



10 20 
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Figure 1. Overlaid polarisation profiles at 1.375 (black lines) and 3.1 GHz (red). In the lower panel of each figure, the solid lines denote 
total power, doted lines linear polarisation and dash-dotted lines circular polarisation. In the upper panel, the absolute PAs are shown 
for 1.375 and 3.1 GHz, following the same colour convention. Also, the error-bars with riser lines (blue) denote the difference between 
the two PAs, where there is a measurement at both frequencies. 
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Figure 1. - continued. 
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Figure 1. - continued. 
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Figure 1. - continued. 
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Figure 2. Same as Figure 1. The top two plots show PSR J1056— 6258, the profile alignment performed on PA (left) and total power 
(right). The bottom two plots show the same for PSR J1359— 6038. The offset between the two alignment methods is 10.7 ms for 
J1056-6258 and 0.7 ms for J1359-6038. 
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at 1.375 than at 3.1 GHz. To further facilitate comparisons, 
we have normalised the flux density at both frequencies with 
respect to the peak of the profile. In doing so, the relative 
spectral dependence of the various components of the pulse 
becomes evident. In most of the pulsars presented here, the 
strongest profile component is the same at both frequencies. 
In these cases, the normalisation also permits direct com- 
parisons between the fractional polarisation of these com- 
ponents. We also draw a comp arison with previous data at 
400 M Hz faamilton et alll977l). 600 MHz jMcCulloch et all 
Il978t) 800 and 950 MHz Ivan Ommen et al.lll997Tl and 1612 
MHz jManchester et alJll98oF where appropriate. 

J0738 4042 (B0736 40). The most striking feature 
of this pulsar is the smooth PA swing, which is broken twice 
by orthogonal jumps at the leading and trailing side of the 
profile at both 1.375 and 3.1 GHz. In fact there is little 
evolution of the polarisation profile between these frequen- 
cies, despite an increase of the total intensity of the leading 
component compared to the main peak of the profile at 3.1 
GHz. The orthogonal jump in the leading component occurs 
exactly at the longitude where both linear and circular po- 
larisation are zero and is accompanied by a change in the 
circular polarisation sense at both frequencies. This is also 
true for the trailing component jump, although the change 
in handedness of the circular polarisation is not as promi- 
nent. Drawing a comparison between our profiles and pre- 
vious profiles at other frequencies, we focus on the leading 
component. The pulse longitude at which the first PA jump 
occurs shifts from the notch between the leading and middle 
component at 631 MHz to a slightly earlier longitude at 950 
MHz. At that frequency, a second local minimum of linear 
polarisation appears at the pulse longitude of the notch. Sur- 
prisingly, the 1612 MHz profile shows the PA in the region 
between the local minima in linear polarisation following the 
previous swing and the jump occurring at the notch, sim- 
ilar to the 631 MHz data. In the 1.375 GHz profile shown 
here, the PA agrees with the 950 MHz data. The most plau- 
sible picture therefore, involves a change in dominant OPM 
between the two local minima of linear polarisation, occur- 
ring between 631 and 950 MHz and subsequently shifting 
the pulse longitude of the PA jump at frequencies above 950 
MHz. 

The comparison of PAs at 1.375 and 3.1 GHz shows that 
not all three segments of the PA profile align. Depending on 
the choice of RM, either the first or the second and third 
parts match. We have chosen the latter in Figure 1, and the 
differences in PA in the second and third segments are zero 
within the uncertainties. Note that the alignment method 
results in alignment in the total power of the main peak 
despite the misalignment of the linear polarisation. 

J0742 2822 (B0740-28). The total power profiles 
at 1.375 and 3.1 GHz have the same overall structure. The 
leading part is made up of two overlapping components, the 
middle part of three and the trailing part of a narrow and 
a broad compon ent (seven components are also suggested 
in iKrameil 1 1994) . However, the ratios of the component 
strengths are different at the two frequencies. In the above 
order, the first, fourth and fifth components have flatter 
spectra than the second and third components. Arguably, 
the same is true for the two trailing components with re- 
spect to components two and three. Both profiles are highly 
linearly polarised with a small amount, right-hand circular 



polarisation across the pulse, in agreement with earlier ob- 
servations. The only obvious difference in the PA profiles is 
a 90° jump in the penultimate component at 3.1 GHz, not 
seen previously at other frequencies. The best alignment in 
PA also aligns the t o tal po wer. 

iGould fc Lvnei (Il998fi have observed this pulsar in 
full polarisation at frequencies up to 1642 MHz, while 
Ivon Hoensbroech fc Xilourisl <ll997bl) have provided profiles 
at 4.85 GHz and 10.55 GHz. The 4.85 GHz profile does not 
demonstrate the orthogonal PA jump in the trailing com- 
ponent and the 10.55 GHz profile lacks a measurement of 
position angle in this component due to the low S/N of 
the linearly polarised flux. Careful investigation of the fre- 
quency dependence of the polarisation in the trailing compo- 
nent will be useful in determining whether there is a mono- 
tonic frequency behaviour of the strength of each OPM. If 
so, there shoul d be a frequency above w hich the dominant 
mode switches jKarastergiou et alJl2005l) . Also, in the trail- 
ing component of the 3.1 GHz profile, the mean circular 
polarisation is very small across the component. Similar be- 
haviour (i.e. a drop in the mean linear polarisation, orthog- 
onal PAs and zero mean circular polarisation) has been ob- 
served in a number of pulsars which show wide distributions 
of circular polarisation, with many instances of high left- 
and high right-hand circular polarisation , resulting in a con- 
siderably higher value for the mean \V \ llKarastergiou et alJ 
2003|; iKarastergiou fc Johnstonll2004l) . We predict that sin- 
gle pulse observations at high frequencies will reveal this 
effect in this pulsar. 

J0835-4510 (B0833 45). The profile of the Vela 
pulsar at 1.375 and 3.1 GHz is made up of three compo- 
nents: two, overlapping at the leading edge and a weaker 
one at the trailing edge. The flux density ratio of these three 
components with respect to each other is different at the 
two different frequencies, in that the second and third com- 
ponents become more prominent at 3.1 GHz. This pulsar is 
very highly linearly polarised at these frequencies, as it is 
at lower frequencies. The PA swing of Vela was the original 
swing used to advocate the RVM. The PA profiles at 1.375 
and 3.1 GHz are virtually identical. Small differences arise 
mainly around the profile centre, where the aforementioned 
changes in relative intensity of the overlapping components 
also occur. 

J0837-4135 (B0835 41). The profile consists of a 
central component and two outriders. The outriders are 
hardly detectable at lower frequencies but become more and 
more prominent at 1.375 and 3.1 GHz. The PA swing in the 
outriders is similar at both frequencies, with a difference of 
zero within the errors. In the middle component, there are 
two peaks of linear polarisation at both frequencies, with 
different ratios. The largest disagreement in PA also occurs 
within this range of the profile. An abrupt PA discontinuity 
at 1.375 GHz can be seen between the leading and mid- 
dle component, despite a much smoother transition at 3.1 
GHz. The trailing peak of linear polarisation in the middle 
component shows a flat, identical PA at both frequencies. 
Also, there is less overall circular polarisation in the mid- 
dle component of the 3.1 GHz profile with respect to 1.375 
GHz. The jump in the 1.375 GHz profile is also seen at lower 
frequencies. 

J0942-5552 (B0940 55). The three main compo- 
nents of this pulsar show the archetypal behaviour of out- 
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riders coming up with respect to the middle component at 
higher frequencies. The PAs at 1.375 and 3.1 GHz are iden- 
tical across the pulse profile. The central component at 3.1 
GHz seems to lag its 1.375 GHz counterpart. The profile 
of linear polarisation in the middle component aligns well, 
despite an extra peak at 3.1 GHz that does not change the 
overall width of the component. 

J1056 6258 (B1054-62). The profiles of this pulsar 
at 1.375 and 3.1 GHz are shown in the top row of Figure 
2. The graph on the left shows the alignment preferred by 
the PA matching routine, which results in an obvious off- 
set between the total intensity profiles. On the right, the 
total intensities have been matched. The offset between the 
two ca ses is 10.7 ms. Using the model of Blaski ewicz et al] 
Jl99lfl . it is possible to convert this delay into a difference in 
an emission height, without obtaining the absolute emission 
height at each fre quency, according to the simple equation 
jMitra fc Lill2004ft 

AR&--xAt (5) 
4 

where c is the speed of light. For this pulsar, we obtain an 
emission height difference of ~ 800 km, which is unlikely. 

Another possible cause of the alignment discrepancy is 
that the overlapping linearly polarised components change 
their relative intensity at 3.1 compared to 1.375 GHz, 
thereby intrinsically changing the PA profile. The change 
is such, that the best PA alignment results in the offset in 
the total power profile. This explanation then suggests that 
the total power profiles aligned between the two frequencies 
is correct and not the PA alignment. The RM computed 
is then different for the two alignment methods; using the 
PA alignment, the RM is —1.0 rad m -2 , where as the total 
power alignment results in an RM of 6.0 rad m -2 , which is 
much closer to the previous published value of 4.0 rad m -2 
JCosta et alJl!99 lft . 

Comparing the profiles of this pulsar, reveals a large re- 
duction in the linear polarisation at 3.1 GHz, especially in 
the leading part of the profile. A 90° jump in PA is seen at 
3.1 GHz lending support to the h ypothesis that the depolar - 
isation is due to competing OPM llKarastergiou et alj200Sf) . 
Also, there is almost no integrated circular polarisation at 
either frequency. 

J1157-6224 (B1154-62). At 3.1 GHz there are 
prominent outriders in the total power profile compared to 
the 1.375 GHz profile. The linear polarisation in the middle 
component is marginally less at 3.1 GHz than at 1.375 GHz. 
The PA swing of the both profiles is complex but identical 
for the pulse longitudes where the linear polarisation is high. 
The circular polarisation shows an S-shaped swing around 
the centre of the pulse, often associated with emission orig- 
inating from close to the magnetic axis jRankinlll986l) . de- 
spite the flat PA profile. At low frequencies this pulsar has a 
single-component profile with some linear polarisation and 
a flat PA profile. 

J1243 6423 (B1240-64). There is more outlying to- 
tal power emission on either side of the central components 
of this pulsar at 3.1 than at 1.375 GHz. A decrease in lin- 
ear polarisation with frequency is also evident, despite both 
profiles having fairly large, right-hand circular polarisation. 
The PA profile is characterised by two discontinuities: the 
first, at pulse longitude 12° appears to be orthogonal at 



1.375 GHz and only ~ 60° at 3.1 GHz; the second, just be- 
fore pulse longitude 20° occurs at a very steep part of the 
PA swing and is not orthogonal at either frequency. The 
main differences in the absolute PAs between 1.375 and 3.1 
GHz are concentrated in a window between pulse longitudes 
16° and 20° which coincide with the steepest parts of both 
PA curves. It is interesting that the general shape of the PA 
profile at both frequencies resembles a v, with a flat segment 
on either side, allowing for 90° discontinuities. In previous 
observations the jump at the trailing edge is seen for the 
first time at 1612 MHz, whereas it is not seen at 950 MHz 
or 631 MHz. 

J1326-5859 (B1323 58). At 1.375 GHz, the profile 
of this pulsar consists of a weak component, followed by 
a strong middle component which trails off, hinting at the 
presence of a weak trailing component. The linear polari- 
sation, which is moderately high across most of the pulse, 
shows a minimum under the leading component, coinciding 
with a 90° discontinuity in the PA. There is a swing in the 
circular polarisation, from right to left hand, in the centre 
of the pulse. At 3.1 GHz, the total power profile looks more 
complicated. The PA profile resembles the lower frequency 
profile, with some notable differences. There is an orthogo- 
nal PA jump at the leading part of the profile, as at 1.375 
GHz. Both PA profiles feature kinks and bends, particularly 
under the middle total power component. The linear polar- 
isation there consists of two overlapping components, more 
clearly seen at 3.1 than at 1.375 GHz and the PA swing is 
also more perturbed at the higher frequency. There is also 
a small right-hand circular polarisation feature just before 
the change in handedness in the 3.1 GHz profile. Also, the 
smooth positive slope in PA joining the middle and trail- 
ing parts of the 1.375 GHz profile becomes what appears to 
be an orthogonal PA transition after a minimum in linear 
polarisation at 3.1 GHz. The total power components have 
different evolution with frequency, the leading component 
having the flattest spectrum. The second component of the 
central pair, which has a steeper spectrum with respect to 
its neighbours, also seems to de-polarise significantly. 

J1327 6222 (B1323 62). The profile of this pulsar 
consists of at least four components, which are more distin- 
guishable at 3.1 than at 1.375 GHz. There is little linear and 
circular polarisation in both profiles. The PA profile at the 
two frequencies does not resemble the S-shaped curve of the 
RVM. In fact, between pulse longitudes 5° and 13° the PA 
profiles are v shaped. At pulse longitude 13° an orthogonal 
jump occurs at both frequencies, whereas at longitude 5° 
the steep increasing PA at 1.375 GHz becomes an apparent 
90° discontinuity at 3.1 GHz. The PAs at 1.375 and 3.1 GHz 
agree over approximately half the profile, namely the lead- 
ing and trailing ends. Significant deviations can be found in 
the middle part of the profile. 

Out of the two strongest components, the leading one 
has a comparatively flatter spectrum and more linear polar- 
isation at 3.1 than at 1.375 GHz. The trailing edge of the 
profile also appears more fractionally polarised at 3.1 than 
at 1.375 GHz. 

J1356 6230 (B1356 60). At 1.375 GHz, the profile 
of this pulsar consists of a strong middle component flanked 
by outriders. Linear polarisation is moderate and a large 
amount of right-hand circular polarisation is seen in the 
middle of the pulse. At 3.1 GHz, the outriders are stronger 
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than the middle component, and the degree of polarisation 
is higher in most of the pulse. The circular polarisation com- 
ponent seen at 1.375 GHz is not present at 3.1 GHz. As with 
J0738— 4042, the PA profile at both frequencies consists of 
three segments, only two of which align between 1.375 and 
3.1 GHz. The PA under each of the three components does 
not vary, with a similar value in the two outer components, 
displaced with respect to the flat PA in the middle compo- 
nent. 

J1359 6038 (B1356 60). This pulsar is shown in 
the second row of Figure 2. Similar to PSR J1056— 6258, the 
total intensity profile at 3.1 GHz lags the 1.375 GHz profile 
when attempting an alignment based on PA, as shown in 
the left plot. To align the total intensity profiles as shown 
in the right plot, a shift of 0.74 ms is required. Equation |S] 
translates this delay into to ~ 55 km of difference in emission 
altitude. However, the additional component coming up at 
3.1 GHz at the trailing edge of the pulse may be responsible 
for the PA differences observed there when aligning the total 
power. This pulsar is highly linearly and circularly polarised, 
with increasing circular polarisation at 3.1 GHz. The PA 
profiles also look identical between the two frequencies. 

J1602-5100 (B1558-50). At 1.375 GHz, this pulsar 
has two main peaks made up of multiple components. The 
PA profile shows two flat segments, joined in the middle by 
a very steep swing, that appears to continuously sweep over 
225°. There is an orthogonal PA jump at 1.375 GHz at the 
very leading edge of the profile (pulse longitude of 12°). At 
3.1 GHz, the last component of the trailing peak becomes 
more prominent in comparison with the other components, 
as does the linear polarisation in the same component. De- 
spite the large complexity of the PA profile at 1.375 and the 
few significant PA measurements at 3.1 GHz, the PAs mostly 
match well, especially at the leading and trailing edges of the 
pulse. 

J1644-4559 (B1641-45). The PA profiles of this 
pulsar resemble the v shaped profiles of PSR J1327— 6222. 
This pulsar is moderately linearly polarised, the fractional 
linear polarisation increasing with frequency. There is an or- 
thogonal PA jump at the leading edge of the profile, seen at 
both 1.375 and 3.1 GHz, whereas an additional 90° PA jump 
is seen at the trailing edge of the 3.1 GHz profile, proceeded 
by a non- instantaneous transition of ~ 90°. The leading and 
trailing parts of the PA profile are flat and do not simulta- 
neously align. In the middle of the profile, although the two 
PA curves have the same overall swing, the kinks are more 
pronounced at 3.1 than at 1.375 GHz. There is also a change 
in the sense of the circular polarisation around pulse longi- 
tude 26° . Previous observations at lower frequencies show 
mainly the effects of scattering on the profile. 

J1752 2806 (B1749-28). At both frequencies, the 
profile of this pulsar shows little linear polarisation and cir- 
cular polarisation that swings from right to left hand. The 
best alignment of PAs between the two frequencies occurs 
at the leading part of the profile, where the PA has a steep 
negative slope. However, the PAs disagree after pulse lon- 
gitude 10° . This coincides with the part of the pulse that 
has a steeper spectral index. An orthogonal PA jump near 
the leading edge of the profile is observed at 631 MHz and 
950 MHz in previous profiles, but not at 400 MHz and 1612 
MHz, although there is very little linear polarisation to make 
a PA measurement. From the 1.375 GHz profile shown here 



we suspect that the steep negative PA slope at the leading 
part of the profile can give the impression of an orthogonal 
jump in observations with poor temporal resolution. 

J1807 0847 (B1804 08). This is another case of a 
profile consisting of a central component with outriders com- 
ing up at 3.1 GHz. Low to moderate linear and circular po- 
larisation are seen across the profiles. The minimum in linear 
polarisation in the leading component at both frequencies, 
occurs at the same time as an orthogonal PA jump. At 1.375 
GHz there is a second orthogonal PA jump in the trailing 
edge of the profile, which can also be inferred from the 3.1 
GHz profile. A systematic deviation in PA between the 2 
frequencies occurs in the middle part of the profile. At 950 
MHz, 90° PA jumps in the leading and trailing edge of the 
profile can be seen, much like the profiles shown here. 

J1825 0935 (B1822-09). As is the case with the 
previous pulsar, there is little change in the polarisa- 
tion profile between 1.375 GHz and 3.1 GHz (other po- 
larisation profiles can be found i n iGould fc Lvnd Il998l . 
von Hoensbroech fc Xilouris! Il997bl) . The leading compo- 
nent remains very highly linearly polarised and the kinks 
in the PA profile of the trailing component coincide with a 
dip in the linear polarisation. 



5 DISCUSSION 

Of the pulsars studied here, only Vela, PSRs J0742-2822, 
J1056-6258 and J1359-6038 show a smooth, unbroken PA 
swing expected under the RVM. RVM fitting can be per- 
formed with confidence only in the former two ca ses. For 
10 of t he pulsars in our list, a has been derived in iRankinl 
( 199330), by using the width of the central component. The 
derived values are mostly in agreement with g eometrical val- 
ues presented bv iLvne fc Manchester! il988l) . who instead 
made use of the RVM. However, as is typical of the pulsar 
population generally (especially at high frequencies), com- 
puting geometrical parameters such as a and (3 must neces- 
sarily be subject to large uncertainties. We therefore do not 
concentrate on the geometrical implications of these results; 
rather we focus on possible causes of PA variations and their 
frequency dependence based purely on a direct comparison 
between the PA curves. We have made note of observations 
at other frequencies but caution that observations made with 
low time resolution can often be misleading. 

We do note however, that many of the pulsars in our 
sample show complex (multi-component) profiles and large 
overall PA swings. This likely indicates that /3 is small in 
these pulsars and the sightline traverses close to the mag- 
netic pole. Also, deviations of the PA from smooth RVM 
curves tend to occur mostly in the centre of profiles either 
because components overlap there or perhaps for some phys- 
ical reason assoicated with core emission. Observations at 
even higher frequencies, where individual components are 
narrower and the core emission much less prominent would 
then reveal smoother PA curves, an effect which is seen, at 
least in some pulsars jvon Hoensbroech et al.lll998l) . 

The 17 pulsars presented in Figures 1 and 2 and de- 
scribed in the previous section demonstrate the diversity 
in emission properties. The number of total intensity com- 
ponents, the degree of polarisation and the swing of the PA 
with pulse longitude vary significantly from pulsar to pulsar. 
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Figure 3. Top panel: The difference in PA between 1.375 and 3.1 GHz, as a function of the position within the width of the pulse, where 
PA measurements are possible. 



Superposing the pulse profiles at 1.375 and 3.1 GHz reveals 
differences in the emission properties between these two fre- 
quencies. These differences, noted in the previous section, 
mostly relate to the ratios of the various component ampli- 
tudes within a profile and the varying degree of linear and 
circular polarisation of individual components at the two 
frequencies. The changing component amplitude ratios are 
a manifestation of the different spectral behaviour of differ- 
ent components, a topic extensively discussed recently by 
iKarastergiou et all i2005t) . 

The pulse widths presented in Table 1 reveal little 
change between 1.375 and 3.1 GHz. This is not surpris- 
ing, considering that both these frequencies are relatively 
high. The flaring of the last open magnetic field lines that 
results in broader profiles at lower frequencies does not ap - 
pear to play a significant role here. iMitra fc Rankinl £2002) 
also point out that components that originate from mag- 
netic field lines closer to the magnetic axis have the same 
widths across wider frequency ranges. It is therefore evi- 
dent that information on the difference in emission heights 
cannot be gleaned from the pulse widths. For the pulsars 
of Figure 2, namely J1056-6258 and J1359-6038, the PA 
alignment method may not be accurate, due to intrinsic 
PA changes. For the remaining pulsars we observed, the PA 
method aligns the profiles to ±1 data bin, the sampling time 



of which is equal to the pulse period over 1024. This places 
a limit on the range of emission height differences between 
1.375 and 3.1 GHz for our pulsars, ranging from ~ 6.7 km for 
the Vela pulsar to ~ 64.8 km for J1602— 5100, as computed 
from Equation |S] 

There are also 4 pulsars in our source list 
(PSRs 0942-5552, 1157-6224, 1356-6230 and 1807-0847) 
that have profiles with an outer component pair flanking 
a central component. In all 4, the relative intensity of the 
outer components with respect to the central component 
is higher at 3.1 than at 1.375 GHz. However, none of 
these sources show spectacular effects when it comes to 
the PA evolution with frequency. In fact, the first three 
pulsars show flat PAs across the central component at 
both frequencies and only PSR J1807— 0847 has a steep PA 
swing in the middle of the profile. Unfortunately, however, 
the linear polarisation is very low in the central component 
of this pulsar, making the frequency comparison close to 
impossible. 

For these 4 pulsars it is possible to derive absolute emis- 
sion heights at each frequency, just by measuring the off- 
set of th e central peak from the middle of the outer com- 
ponents iGangadhara fc GuptallioOltlDvks et aill2004h . Al- 
though the direction of the offset is in agreement with the 
fact that the central component is emitted at a lower height 
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(the leading outer component is always closer to the central 
component than its trailing counterpart), it is clear that the 
offset is identical at the two frequencies for all 4 of the afore- 
mentioned pulsars. This constitutes further evidence that 
the emission heights are very similar between the two ob- 
served frequencies. The absolute emiss ion heights derived, 
using equation 7 of iDvks et all l|2004f) . are 14.5, 80, 33.5 
and 4.5 km respectively for the aforementioned 4 pulsars. 
PSR 1807—6230 has the shortest rotation period from this 
selection, limiting the radius of its magnetosphere. 

In the PA profiles, there is an obvious overall similarity 
between the two frequencies. Occasionally, the PA is differ- 
ent in narrow or extended parts of the profile, although even 
then, the agreement between frequencies outweighs the dis- 
agreement. To illustrate this, we have summarised in Figure 
3 the PA differences dPA in all the pulsars we observed. 
We initially identified the pulse widths over which PA mea- 
surements were possible, for all the pulsars. The abscissae 
on the top panel denote the pulse longitude as a fraction 
of this width. Note that these widths are sometimes greater 
than the Wio widths of Table 1. This creates a normalised 
pulse longitude scale for all the pulsars. The points on the 
plot represent the difference in PA, or dPA, as shown in the 
plots of Figures 1 and 2, plotted against their location within 
the pulse. Furthermore, the second panel in Figure 3 shows 
the histograms of dPA in each of three regions of the top 
panel. The intention is to demonstrate that the distribution 
of dPA does not obviously depend on the location within 
the pulse. The standard deviation of all the dPA values in 
the top panel is 16.6°, compared to 13.9°, 13.8° and 21.2° in 
the three regions of the bottom panel. The higher standard 
deviation value of the third region is partially caused by the 
large negative dPA values in PSR J1644— 4559. The average 
error in dPA is ~ 7, so the standard deviations measured 
are between 2 and 3 times this value, which consolidates the 
observation that the PA differences cannot merely be put 
down to noise. 

The histograms in the bottom row of Figure 3 can be 
used to identify OPM as a possible source of PA disagree- 
ment between 1.375 and 3.1 GHz. If at some intermediate 
frequency, the dominant mode of polarisation has changed, 
the PA at the high frequency will be orthogonal to the low 
frequency. Interestingly, in Figure 3 there appears to be a 
Gaussian spread of dPA around its mean of zero. The in- 
stances of dPA close to 90° or —90° are few and do not rep- 
resent a clear sub-group of the dPA distributions. This fact 
suggests that for our sample of sources, changes in the dom- 
inant OPM are not common between these two frequencies. 
The general decrease in the linear polarisation, however, is 
an indicatio n that the strengths of th e OPMs are becoming 
more equal jKarastergiou et al.ll2002l) . 

The measured standard deviation in the distribution of 
dPA places constraints on the frequency dependence of pos- 
sible mechanisms that set the PA. If the superposed modes 
of polarisation are not orthogonal and have different spectral 
behaviour, a dependence of the PA on frequency will arise. 
The different spectral behaviour is evident by the change in 
linear polarisation between the two frequencies. The left plot 
in Figure 4 shows the effect of superposed, non-orthogonal 
modes of polarisation on the observed degree of linear po- 
larisation and PA. From left to right, the lines correspond 
to a deviation from orthogonality of 1°, 3°, 5°, 10° and 20°. 



The vertical line has a length of 13°, as a measure of the 
standard deviation of dPA. This figure demonstrates how 
large dPAs combined with small changes in the degree of 
linear polarisation, require larger departures from orthogo- 
nality of the polarisation modes. It also demonstrates that 
non-orthogonality of the modes results in the minimum de- 
gree of linear polarisation not being zero, as would be the 
case with orthogonal modes. More importantly however, this 
figure shows that in cases where the degree of linear polari- 
sation changes significantly from one to the other frequency, 
non-orthogonality will lead to larger values of dPA. We con- 
ducted the following test: for all bins where the fractional 
linear polarisation is greater than 0.5 at 1.375 GHz, we mea- 
sure the difference of the degree of linear polarisation and 
dPA between the two frequencies. The results are shown in 
the right plot of Figure 4. There is a hint that the points 
which correspond to a large change in the fractional linear 
polarisation, also correspond to larger dPA. The correlation 
however is very weak, which demonstrates that an explana- 
tion involving non-orthogonal polarisation modes is unlikely. 

In the individual plots of Figures 1 and 2, we identify 
two types of deviation of dPA from zero. The first consists of 
a part of the profile where the PAs at the two different fre- 
quencies deviate in a quasi-random way, resulting in a vari- 
able dPA which gradually returns to zero. Good examples 
of this behaviour are PSRs J1243-6423 and J1327-6222. 
In the second type dPA has a constant but non-zero value 
across a fraction of the pulse profile. The PAs at 1.375 and 
3.1 GHz run parallel to each other and then converge to 
being identical. Examples of this can be found in PSRs 
J0738-4042 and J1644-4559. 

In the pulsars that belong to the first type, the total 
power profiles are complex with generally low linear polar- 
isation. We can often discern what appear to be several, 
partially overlapping components of linear polarisation, co- 
inciding with the most perturbed parts of the PA swing. 
The question then arises whether the observed changes in 
PA are merely due to changes in the relative strengths of 
overlapping components of linear polarisation. It seems ap- 
propriate to conduct such tests for PSRs J1243— 6423 and 
J1327— 6222, to formulate a model. From our data, the ratios 
of the various components in e.g. PSR J1243— 6423 are cer- 
tainly different at 1.375 and 3.1 GHz, however single pulse 
data is required to properly acc ount for the effects of su - 
perposed modes of polarisation jKarastergiou et al.l [2003). 
Also, in the case of the Vela pulsar where the linear polar- 
isation is very high, the second component in the leading 
edge of the profile becomes more prominent at 3.1 GHz and 
at the same time the PA profile is distorted at the same pulse 
longitudes. An interpretation based on overlapping compo- 
nents with intrinsically different PAs again seems likely. 

For the second type of dPA offsets, the same reasoning 
cannot explain the observations. It is tempting to attribute 
the fact that some components show a constant offset in PA 
between 1.375 and 3.1 GHz, e.g. the leading component of 
PSR J0738— 4042, to a different Faraday Rotation measure. 
Although we cannot rule this possibility out, there seems to 
be little support in terms of the underlying physics. How- 
ever, a proper test requires one more observing frequency, 
to check if the dependence of PA on frequency is of the same 
nature as Faraday Rotation. The fact that the observed dif- 
ference in PA corresponds only to a few units of rotation 
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Figure 4. Left: The effect of superposed, non-orthogonal polarisation modes on the observed PA and degree of linear polarisation (L/I). 
The different curves, from left to right, correspond to a deviation from orthogonality of 1°, 3°, 5°, 10° and 20°. Right: |dPA| plotted 
against the difference in the fractional linear polarisation, for all instances of high linear polarisation at 1.375 GHz. 



measure does not permit us to conduct this test within one 
of the bands of observation. 



6 CONCLUSIONS 

We have presented for the first time superposed polarimetric 
profiles of 17 pulsars at 1.375 and 3.1 GHz, with absolute 
values of PA. The presentation demonstrates the frequency 
dependence of the PA as a function of pulse longitude. We 
find that, in general: 

• The PA profiles are similar at the two frequencies, de- 
spite the lack of resemblance to the RVM curve. 

• We can identify a tentative association between the 
changes in the PA between 1.375 and 3.1 GHz and the 
change in relative amplitude of overlapping total power com- 
ponents. This observation clashes with the RVM in that the 
polarization of the emission at each pulse longitude is con- 
strained by the pulse component in which it belongs, rather 
than being determined purely by geometry. 

• The similarities permit an alignment method based on 
the PA profile, which allows us to calculate very accurate 
RMs. The method performs very well for all but two pulsars: 
J1056— 6258 and J1359— 6038, where an explanation related 
to intrinsic development of the PA profile seems more likely. 

• In our sample, there are not many instances of the PA 
at one frequency being orthogonal to the other. 

• Using our measurements of the fractional linear po- 
larisation, we place constraints on the degree of non- 
orthogonality between the polarisation modes. 

• Although in most pulsars the differences in PA between 
the two frequencies vary on short scales of pulse longitude, in 
some pulsars such as J0738— 4042, we have found a constant 
offset in PA across extended parts of the profile. We intend 
to perform specific observations on this pulsar to identify 
whether this offset can be accounted for by a different RM 
value, or some other physical effect. 
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